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6 ■ ABSTRACT 

& ■ 

r/3 , We report the detection of the slow moving wind into which the compact 

CG ■ 

supernova remnant SN 1997ab is expanding. Echelle spectroscopy provides 
clear evidence for a wellresolved narrow (Full Width at Zero Intensity, FWZI ~ 
?—i ' 180 kms -1 ) P-Cygni profile, both in Ha and H/3 , superimposed on the broad 

emission lines of this compact supernova remnant. From theoretical arguments 
we know that the broad and strong emission lines imply a circumstellar density 
(n > 10 7 cm' 3 ). This, together with our detection, implies a massive and slow 
stellar wind experienced by the progenitor star shortly prior to the explosion. 

Key words: Supernovae and Supernova Remnants: general - Supernovae 
and Supernova Remnants: SN 1997ab - Circumstellar Medium. 



* Based on observations made with the 4.2-m WHT operated on the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. 
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1 INTRODUCTION 

SN1997ab, like SN1987F or SN1988Z (Filippenko 1989, Turatto et al. 1993) is a peculiar 
type II supernova (SN) with a luminosity of Mb = —17.5 one year after its discovery (Hagen, 
Engels & Reimers 1997). It was first detected on April 11, 1996, in the dwarf galaxy HS 
0948+2018, far from the galaxy center, showing broad Balmer emission lines, FWHM ~ 2500 
to 3000 kms -1 , without the typical P-Cygni profiles of type II SN. The redshift derived 
from the narrow emission lines seen in SN 1997ab is 0.012 (Hagen et al. 1997). These peculiar 
SN, also called "Seyfert 1 impostors" (Filippenko 1989), given their profound resemblance to 
the emission spectra from Seyfert 1 nuclei, are sometimes also identified as radio supernovae 
(RSN; see Van Dyk et al. 1993, 1996). The standard explanation for these sources (Chugai 
1990, Terlevich et al. 1992, 1995) demands the expansion of a supernova shock into a dense 
circumstellar medium, and thus it is believed that what we see is the SN remnant and not the 
SN itself. Indeed, the release of 10 51 erg into a high density medium (> 10 6 cm -3 ) is known 
to lead to a rapidly evolving remnant, which skips its quasi-adiabatic Sedov evolutionary 
phase to become strongly radiative, reaching luminosities > 10 9 L Q even before the ejecta 
is fully thermalized at the reverse shock (see Wheeler et al. 1980; Terlevich et al. 1992). 
The rapid evolution is a direct consequence of the high densities, which promote the earlier 
onset of strong radiative cooling in the matter swept up and strongly accelerated by the 
blast wave. Strong cooling causes the collapse of the swept up gas into a thin shell while 
still moving at several thousands of kms -1 and with it, to an outburst of radiation capable 
of ionizing the cooling shell, the unshocked ejecta and the background gas. In this way the 
supernova energy, the 10 51 erg originally released as kinetic energy, is radiated away in just 
a matter of a few years. This causes the most luminous SN remnants, although they only 
acquire sizes ~ a few 10 16 cm, and hence have been termed "compact" supernova remnants 
(CSNRs). 

Here we provide unambiguous evidence for the existence of a slowly expanding, dense 
circumstellar medium around SN 1997ab and work out its relevant physical properties. 

2 THE OBSERVATIONS 

On May 30 1997, the supernova remnant, SN1997ab (Hagen et al. 1997), was observed with 
the Utrecht Echelle spectrograph at the 4.2-m William Herschel Telescope at the "Roque de 
los Muchachos" observatory (La Palma, Spain). A 1200s exposure of the supernova remnant 
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was obtained using a 2148x2148 Tek CCD and the 79.0 lines/mm echelle, reaching a S/N 
~ 30 at Ha. The orders covered are 45 to 30, from ~ 4130 A- 7340 A, including H/3 and 
Ha . The spectral resolution of our setup was 0.13 A in the H/3 region and 0.18 A in the 
Ha region. This led to a resolution of ~ 8 kms -1 . The CCD frames were reduced using 
standard procedures with the tasks in IRAF. The data was debiased, trimmed and flat- 
field using a normalized flat field produced with apflatten. The spectra were extracted using 
doecslit with appropriate parameters. Wavelength calibration was performed by using the 
comparison spectrum of Thorium- Argon arcs. Flux calibration was accomplished through 
the observation of standard stars. The atmospheric extinction correction was also applied at 
the time of calibration using the mean extinction curve for La Palma. No redshift correction 
was applied, and therefore the wavelengths are the observed ones, unless otherwise stated. 
The seeing was less than 1.5 arc seconds. 

2.1 Line analysis 

Figure [I] shows the spectral regions around Ha and B./3 lines. Both of these lines have two 
components: a broad component with a FWHM ~ 1800 kms -1 and a narrow (FWZI ~ 
180 kms -1 ) P-Cygni profile superimposed on the broad lines. Note that the broad compo- 
nents present a flat-topped profile, typical of expanding shells, and that the P-Cygni narrow 
lines are somewhat displaced towards longer wavelengths. The broad lines indicate a large 
expansion velocity of (at least) ~ 2000 to 4000 kms -1 , as measured from their red and 
blue wings, respectively. However, the fact that the P-Cygni profiles are shifted implies that 
the red part of the broad components of both, Ha and H/3 , are probably self-absorbed and 
therefore missing. As a consequence, the intrinsic true width of the broad profiles could be 
almost twice as large as what one can at first glance measure. Also, given the small wave- 
length range covered in each portion of our spectra, the broad Ha and B./3 lines are not fully 
sampled and their wings and adjacent continuum are missing. Therefore, regardless of their 
asymmetry, the FWZI cannot be measured in our Echelle spectra. A better determination 
of the FWZI is to be found in previous published (low-dispersion) spectra of this CSNR 
(Hagen et al. 1997, Salamanca et al. 1998) which indicates a velocity of ~ 6600 kms -1 . 

The narrow P-Cygni lines shown in Figure |I]on top of the broad Ha and B./3 , must arise 
in a medium with a much smaller velocity, v w ~ 90 kms -1 , as deduced from their half 
width at zero intensity. 
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Figure 1. The H/3 and Ha profiles of SN1997ab. The flux is in units of erg A 1 s 1 cm 2 . A "smooth" has been applied to 
both spectra to reduce the noise in the wings. 
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3 DISCUSSION 

The detection of narrow P-Cygni profiles in Ha and H/5 , superimposed on the broad emis- 
sion lines produced by the compact supernova remnant SN 1997ab, proves the existence of a 
slowly expanding, dense medium into which the CSNR shock evolves. From the analysis of 
the broad lines one can derive the physical conditions of the SN ejecta and its main shock. 
On the other hand, the narrow P-Cygni lines tell us about the density and extent of the 
CSM. One word of warning however: some of the narrow Balmer emission is contaminated by 
the extended emission from the host galaxy, a starburst galaxy of low excitation. Neverthe- 
less, the emission coming from the host galaxy is very weak, and therefore its contamination 
is negligible in a first approximation. More detailed analysis will be referred to uncoming 
publications (Salamanca et al. 1998). 

To estimate the density and extent of the envelope, we assumed that the wind from the 
progenitor star that created it is freely expanding with a velocity, v w . The density (p) of the 
slow wind then falls as r -2 , under the assumption of a constant mass-loss rate, M, 



3.1 The broad lines and the CSNR shock 

When the shock in the circumstellar medium is the dominant source of luminosity, rather 
than the release and escape of stored radioactive energy, then the luminosity in the broad 
Ha line is simply proportional to the kinetic energy dissipated per unit time across the shock 
front (Terlevich 1994, Chugai 1991), L s = 1/2M C smv 2 , where v s is the shock velocity and 
Mcsm — —Vs is the rate at which the circumstellar medium (CSM) is swept up. The r -2 
density distribution will ensure that the shock velocity can be regarded as constant (until 
an appreciable fraction of the kinetic energy is radiated away). This leads to 



where we have taken into account the fact that only half of the radiation from the gas 
cooling behind the shock will contribute to the ionization of the shell, while the other half 
will be transported outwards. The value of the efficiency factor sua, depends on the time 
elapsed since the explosion of the supernova, and ranges from a maximum value ~ 0.1 when 



M = 4nr 2 v w p 



(1) 



t Broad _ 1 ^ 3 

4 v w 



(2) 
(3) 
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catastrophic cooling starts, and then steadily drops to nearly zero within the evolutionary 
time of the CSNR (Terlevich 1994 and Cid-Fernandes & Terlevich 1994). 

In the case of SN 1997ab, v s and v w are easy to determine from the FWZI of the broad 
and narrow lines, respectively: v s is about ~ 6600 kms" 1 and for v w we shall adopt a 
value of 90 kms -1 (see previous section). The luminosity in Ha is ~ 5 x 10 41 erg s -1 . This 
value, taken from Hagen et al. (1997), is probably a lower limit, since we did not take into 
account reddening or obscuration effects. The value of eua is directly related to the Balmer 
decrement, since it depends essentially on the density of the CSM. From previous spectra 
of SN 1997ab (Hagen et al. 1997) we have ~ 6 which led us to the conclusion that e# a 
must be ~ 0.1 (see also Terlevich 1994). Note also that, as SN 1997ab was not seen on April 
5th 1995 (Hagen et al. 1997), then an upper limit to the time elapsed since the explosion up 
to the end of May 1997, is 783 days, which is again consistent with large values of e# a . 

The values of v s , v w , L# a , etc. have been obtained from spectra taken at epochs separated 
by 1 to 3 months. However, we regard this as a valid approximation as these objects fade 
very slowly. 

Substituting these values into equation |2] we obtain a mass-loss rate M(for eu a = 0.1) 
~ 10~ 2 MQyr -1 , which in what follows is used as reference value. 

The radius of the shock, or inner radius (Ri) of the unshocked CSM able to produce the 
narrow P-Cygni profile, can be expressed as a function of the local density. From equations |l| 
and ^| we have: 



Ri = 1.54 x 10 



16 



,1/2 



1.54 x 10 16 



(M /10" 2 MQyr- 1 ) 
{vyj/dOkms- 1 ) 

(L Ha /5 x IQ^ergs- 1 ) 
(rii I 10 7 cm -3 ) 



cm (4) 



10 7 cm 3 

1/2 / «. \ " 3 /2 



6.6 x 10 8 cms~ 



; cm (5) 

where n« is the density of H atoms at Ri and p = IAuhttih (accounting for a 10% He 
abundance by number). 

3.2 Narrow P-Cygni lines and the CSM 

We use the luminosities of the emission narrow Balmer lines to estimate the unshocked mass 
and outer radius of the CSM. These are given by: 

LZfrZ = r ^r 2 n 2 a B hu dr (6) 
JRi 

were ag is the case B recombination coefficient and we have supposed that the CSM emits 
isotropically and that it is fully ionized. 
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For this estimation we use either the luminosity of the narrow Ha or H/3 lines: The 
Ha line is more intense and it therefore has a better signal-to-noise ratio, but in many cases 
it may be enhanced by collisions and therefore, the above relation (|6]) would no longer be 
valid. The H/3 line is not affected by collisional effects, but it is by reddening. The ratio 
Ha /H/3 is ~ 20, which suggests a high density medium, n ~ 10 s - 10 10 , for an ionization 
parameter logU = -2 and -4 respectively. In what follows we will use the observed values of 
the luminosities in Ha and H/3 narrow emission lines. They will gives us an upper and lower 
limit respectively to the density and extent of the CSM. 

The integration of the above equation using n = rii(r / Ri)~ 2 , yields 

• \ 3/2 

a B hv Q ( M\ 1/2 / Ri 



j narrow _ ^ £>, ""° „H * I 1 "M (7\ 

^Balmer ~ 2 ^/\lAm H ) 3 ^ \v w ) * \ R ' UJ 



3.05 x l0 37 a B v o 



(M/lQ- 2 M Q yr- 1 ) 



-i 3/2 



a, 



V 2 / R 



10 7 cm- 3 J V R, 



erg s 1 (8) 



Using the relationship between M and L^ r ° (see equation one gets: 



T narrow Aqt B hv ( Lff° ad \ 3/2 1/2 / _ Rj_\ 

Balmer " nV*(lAm H e Ha )W { v 3 ) n< V R J 1 ' 



3.05 x 10 37 a B u o 



(Lf° ad /5 x lO^er^s- 1 ) 1372 



71: 



l ' 2 (. Ri 



10 7 cm" 3 J \ R 



1 - Tr) erff[&V) 



(vs/QQOOkms- 1 ) 3 

Substituting the appropriate values for the £™£ row or L n ^j 3 row an d aB and v Q we estimated 
the value of the outer CSM radius, R Q , as a function of the density (rij) at the inner edge 
(Ri) of the CSM. The observed value of L n H a ^ row is 4 x 10 39 ergs _1 and of L™ row 2 x 10 38 
erg s^ 1 (Hq = 75 kms~ 1 Mpc~ 1 ). The true value must be certainly larger, since the blue 
side of the line is in absorption. To account for this and possible obscuration effects we have 
multiplied the observed emission luminosities by a factor of 2. The values of a B taken are: 
a Ha = 2.59 x 10~ 13 cm" 3 s" 1 and a H p = 3.07 x 10" 14 cm" 3 s" 1 (Osterbrock 1989). 

To explain the observed narrow H/3 emission, the CSM must have a minimum mass, and 
this imposes a natural lower bound on rif. 



7ii>3x 10 5 cm" 3 ' Hp 



( L«™ \ 2 \(M/lO- 2 M & yr- 1 ) 



\10 38 erg s 1 



(11) 



and thus, a less dense CSM could not produce the observed B.j3 narrow line. This relationship 
is valid for an infinite CSM. For a finite medium, the density must be even larger. 

Tables 1 and 2 lists the values of R iy R Q , the density (n c ) at R Q , as well as the mass of 
the swept up CSM (up to r = Ri) and of unshocked CSM mass (from Ri to R Q ) for different 
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Table 1. CSM values inferred from the H™ arrotu line 
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Table 2. CSM values inferred from the Wl arrow line 
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values of rij. We as well list the values of the duration of the progenitor wind, t W i n d = R D 
/ 90 kms" 1 , and the time for the supernova shock to reach the inner radius, tsN = R, / 
10000 kms -1 , if assumed to travel at 10000 kms -1 . 

Note that the results for both lines agree only for high densities (~ 10 8 cm _3 or bigger), 
as we would expect from the high value of (see previous section). For such high densities, 
the extent of the CSM becomes extremely small, and the density decreases less than one 
order of magnitude only. This would imply also that the wind that created such CSM was 
of very short duration and that the supernova shock reached it on very short times scales 
(less than two months). However, in the derivation of all these parameters, we have assumed 
that the narrow Balmer emission lines are due only to recombination effects, and we did not 
take into account any collisional effects. With the available data we cannot constrain any 
further these parameters. We need more observations to monitor the narrow P-Cygni profile, 
since its disappearance will indicate the end of the dense CSM. The numbers presented in 
tables 1 and 2 must be regarded therefore as an indication of the possible range of values. 
Note finally that SN 1997ab was not detected in X-rays by ROSAT (A. Fabian, private 
communication), in agreement with the presence of a very large column density of neutral 
material surrounding the SN. and thus a massive wind. 

The rapidly evolving supernova remnant is a peculiar emitting source, in the sense that 
all of its radiation arises from matter at very high speeds (several thousands of kms -1 ). 
This, irrespectively of how dense the circumstellar medium may be, can only be absorbed by 
a medium between the source and the observer moving with a similar velocity. On the other 
hand, the P-Cygni profile results from the emission caused by the precursor radiation ahead 
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of the blast wave, combined with the absorption due to screening of the radiative supernova 
remnant shell and the photoionized freely expanding ejecta. From this it follows that a sector 
of the circumstellar medium, a broad rim, overlapping with the section of the high velocity 
shell that expands near the plane of the sky and/or with the photoionized freely expanding 
ejecta also moving near the plane of the sky, is the sector that can absorb it and provoke 
the P-Cygni profiles here shown. In this respect, the P-Cygni profiles of SN 1997ab are very 
different to those caused by a wind in front of a standing source (e.g. Rublev 1960, Kuan & 
Kuhi 1975). A complete fit to the emission profiles is now in preparation (Cid-Fernandes et 
al. 1998). 

The fact that both the emission and absorption P-Cygni profiles are so similar confirms 
also the fact that the absorbing CSM rim is in fact very broad, and it probably spans from 
the edge of the expanding cooling shell to the photoionized inner most section of the freely 
expanding ejecta. In this way also, there is a central neutral ejecta region that must cause 
the absorption of the redest sections of the broad emission lines. 

To our knowledge, a narrow P-Cygni profile atop a broad emission line such as that of 
SN 1997ab has never been detected in any supernova remnant. f\ Broader P-Cygni profiles 
(but still narrower than those seen in typical type II SN) have been detected in other peculiar 
SN, such as SN 1979C (Panagia et al. 1980), SN 1984E (Dopita et al. 1984), SN 1994aj 
(Benetti et al. 1998) or SN 1996L (Benetti et al. 1996). However, the CSM velocities inferred 
from the widths of all of these P-Cygni profiles are ~ 1000 km s _1 . If a similar interpretation 
is given to these lines, then the large values of v s lead for example in the case of 1994aj to 
an M ~ 3.3 10 -3 M m yr _1 and to an Rj = 1.63 x 10 15 ^=i — , i.e. ~ one order smaller than 
the values deduced for 1997ab. Our result thus broadens the range of values to be expected 
for the massive winds that occur prior to the explosion. 
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